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Abstract 

Quarkonium production has been considered as a tool to study the medium formed 
in high energy nuclear collisions, assuming that the formation of a hot and dense 
environment modifies the production pattern observed in elementary collisions. The 
basic features measured there are the relative fractions of hidden to open heavy 
flavor and the relative fractions of the different hidden heavy flavor states. Hence 
the essential question is if and how these quantities are modified in nuclear collisions. 
We show how the relevant data must be calibrated, i.e., what reference has to be 
used, in order to determine this in a model-independent way. 



The original suggestion for charmonium production as a means to test the formation 
of a deconfined medium in high energy nuclear collisions was based on the idea that 
color screening in such a medium would prevent the binding of charm quarks to a 
color neutral J/ip [TJ. The J/ip production process in elementary hadronic collisions 
(taking pp as example) begins with the formation of a cc pair; this pair can then 
either lead to open charm production (about 90 %) or subsequently bind to form a 
charmonium state (about 10 % for all charmonia). A schematic illutration (Fig. [Q 
shows the dominant high energy reaction through gluon fusion. 
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Figure 1: Schematic view of J/ip production in pp collisions 
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The initial cc production can be calculated in terms of the parton distribution func- 
tions f p of the relevant hadrons and the pertubative partonic cross section. The full 
description of charmonium binding has so far resisted various theoretical attempts; 
on the other hand, the process is in good approximation independent of the incident 
hadronic collision energy |2][3]. This is a consequence of the fact that the heavy 
quark propagator in the reaction gg — > cc strongly dampens the mass variation of 
the cc pair with incident energy. Thus the fractions of the produced cc system into 
hidden vs. open charm as well as those for the different charmonium states are ap- 
proximately constant; once determined at one energy, they remain the same also for 
different energies. As a result, the phenomeno logical color evaporation model [4-7J 
provides a good description of charmonium production through the form 



and correspondingly for the other charmonium states. Here the constant g C c^>j/^ 
specifies what fraction of the total cc production cross section goes into J/ip pro- 
duction; in pp collisions it is typically about 2 %. The set of the different constants 
g C c->i for the different charmonium states i thus effectively characterizes charmonium 
production in the absence of a medium. 

A further important aspect of quarkonium production in elementary collisions is 
that the observed (IS) ground states J/ip and T are in both cases partially produced 
through feed-down from higher excited states [SHU]- Of the observed J/ip rates, only 
some 60 % is a directly produced J/ip{lS) state; about 30 % comes from x c (lP) 
and 10 % from ip'{2S) decay. Because of the narrow width of the excited states, 
their decay occurs well outside any interaction region. 

The features we have here summarized for charm and charmonium production are 
readily extended to that of bottom and bottomonium. To simplify the discussion, we 
shall continue referring to the charmonium case, keeping in mind that all arguments 
apply as well to bottomonia. Given the patterns observed in elementary collisions, 
we want to see how they are modified in the presence of a medium, as provided by 
nuclear collisions. From the point of view of production dynamics, one way such 
modifications can arise is as initial state effects, which take place before the cc pair 
is produced. The main possibilities considered so far are nuclear modifications of the 
parton distribution functions (shadowing or antishadowing) and a possible energy 
loss of the partons passing through the nuclear medium to produce the cc. Once 
produced, the pair can encounter final state effects, either in the form of a phase 
space shift already of the cc, e.g., through an energy loss of the unbound charm 
quarks, or through effects on the nascent or fully formed charmonium state. Such 
effects may arise from the passage through the cold nuclear medium, or because of 
the presence of the medium newly produced in the nuclear collision. The latter is 
evidently what we have in mind when we want to use quarkonia to study quark- 
gluon plasma production. The difficulties encountered over the past years in 
arriving at a conclusive analysis of the relevant nuclear collision data on J/ip pro- 
duction are largely due to the problem of parametrizing the different effects and 
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then constructing a convincing model correctly incorporating all of them. We want 
to show here that today experimental means have become available which allow us 
to calibrate the measured results in a way which avoids these difficulties. 

Before we turn to the problem of determining medium effects on quarkonium pro- 
duction, we recall the two main conceptual approaches on what may happen; the 
two can perhaps best be labelled suppression and enhancement. Color screening in 
a quark-gluon plasma will decrease the quarkonium binding, both in strength and in 
its spatial range, and this should for sufficiently energetic nucleus-nucleus collisions 
lead to quarkonium dissociation or melting. Since the larger and less tightly bound 
states will melt at lower temperature or energy density than the ground states, color 
screening will produce sequential suppression [T1|IT5]. We illustrate this for the J /if). 
After an initial threshold melting the ip' and hence removing its feed-down compo- 
nent for J /if) production, there will be a second threshold for Xc melting and then 
finally a third, at which the direct J/ip(ls) is dissociated. The resulting pattern is 
illustrated in Fig. [2j We have here introduced something denoted as J j ip survival 
probability. Theoretically, this is the chance of a J/ip to persist as a bound state in 
a deconfined medium. How to properly define this quantity as a useful observable 
in experimental studies is the subject of this paper. 
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Figure 2: Schematic view of sequential J/ip suppression in a deconfined medium 

The other alternative, J/ip enhancement, assumes, in accord with color screening 
and the resulting suppression, that at sufficiently high energies there is an almost 
complete dissociation of the J/^'s produced in primary nucleon-nucleon interactions. 
On the other hand, at such collision energies, these interactions lead to abundant 
cc production; the rate for this process grows faster than that for the production of 
light quarks, and if the cc pairs remain present in the evolution of the medium, the 
system will at the hadronisation point show an oversaturation of charm, compared 
to the predicted thermal abundance. If these charm quarks have become part of an 
equilibrated medium and as such undergo hadronisation in the form of statistical 
combination, then such secondary charmonium formation can convert more cc pairs 
into J/ip's than the dynamical primary production mechanism, thus leading to an 
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effective J/ip enhancement [TEH!!?] , as shown in Fig. [3J We note that this scenario 
invokes a new and so far unknown binding dynamics to form charmonia; the outcome 
is assumed to be determined simply by the relative abundance of charm quarks [20J . 




Energy Density 



Figure 3: Schematic view of secondary production through statistical combination 

We now come to our central question: how to calibrate the J/ip survival? Since 
we are interested in using quarkonium production as a tool to study the medium 
produced in nuclear collisions, our primary concern is not if such collisions produce 
more or fewer cc pairs than proton-proton collisions, but rather if the presence 
of the medium modifies the fraction of produced cc pairs going into charmonium 
formation. In other words, the crucial quantitity is the amount of charmonium 
production relative to that of open charm [21] . To illustrate: in pp collisions, about 
2 % of the total cc production goes into J / ip's. If in high energy nuclear collisions the 
total cc production rate were reduced by a factor two, but we still have 2 % of these 
going into J/ip's, then evidently AA collisions do not modify J/ip binding. This is, 
of course, strictly true only for fully integrated rates, which is difficult to check; we 
return to the issue of phase space variations shortly. But we note here already that 
in the example just given, the number of J/^'s produced in AA collisions would be 
half of that obtained by scaling the results from pp interactions. That is, however, 
not an indication of any J/ip suppression; it is just the consequence of having fewer 
cc pairs to start with - the formation probability of charmonia has remained the 
same as in pp collisions, 2 %. It is thus crucial to correctly calibrate the survival in 
nuclear collisions. 

To achieve that, we recall: both sequential suppression and statistical enhancement 
state that the total number of cc pairs produced in nuclear collisions is distributed 
among hidden and open charm differently than it is in proton-proton collisions. The 
essential observable is thus the ratio of charmonium states to open charm in nuclear 
collisions, compared to that in pp interactions. Sequential suppression predicts it 
to decrease with centrality and eventually vanish, statistical enhancement has it 
increasing with centrality, in many studies beyond the pp value. Hence the relevant 
observable is the fraction of charmonia to open charm, or more generally, that of 
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quarkonia to the relevant open heavy flavor production [2~TII2"2"] . In this quantity, 
if measured over the entire phase space, the effects of possible initial state nuclear 
modifications - shadowing/antishadowing, parton energy loss - cancel out, so that 
whatever changes it shows relative to the pp pattern is due to final state effects. 

It should be noted here that in actual applications, it will most likely not be neces- 
sary (nor generally possible) to measure the total open charm production rate. The 
relative abundances of the different open charm states produced in high energy colli- 
sions are so far found to be in good agreement with the predictions of the statistical 
hadronization model [221121], with a universal hadronisation temperature of some 
150 - 170 MeV. For open charm production, this model assumes an initial dynamical 
(perturbative) cc production; subsequently, at the hadronisation point, the charm 
quarks thus produced form open charm hadrons according to their statistical abun- 
dances. The relative production rates are thus totally determined in terms of the 
mass of the open charm state and the universal hadronisation temperature. Such 
a description is found to hold very well for open charm production in high energy 
interactions from e + e~ annihilation through pp and pA collisions; studies of the AA 
behavior are underway. However, the model does not work for charmonium produc- 
tion, neither in elementary (e + e~, pp) nor in nuclear interactions, since the binding 
mechanism forming charmonia is evidently of dynamical and not of statistical origin. 
In nuclear collisions, there are in addition medium effects on the binding. 

As a consequence, we expect that the relative abundances of the different open 
charm states will not depend on the collision energy and remain essentially the 
same for elementary and nuclear collisions. Analogous to eq. ([TJ, the production 
rate N(hh — > Di) for a specific open charm state Di is then a constant fraction of 
the total open charm rate N(hh — > cc), 

N(hh -> Di) ~ g c - c ^ Di N{hh -> cc). (2) 

The crucial observable, the rate of (J/-0)/open charm, can thus in good approxi- 
mation be taken as the rate of (J/tp)/Di for a specific D-meson state. This will, of 
course, greatly facilitate the analyis. 

The final state effect of interest to us is that caused by the newly produced medium. 
We therefore have to check to what extent the nascent or fully formed J/ip is already 
dissociated by the cold nuclear matter of the target or the projectile. This can be 
done by studying the ratio hidden to open charm in high energy pA collisions. It is in 
fact known that in such interactions there exist reductions of J/ip production beyond 
the scaled pp results, but the origin of these is not unambiguously clarified. They 
could arise largely from initial state effects; but if there is a stronger suppression of 
ip' than of J/ip at energies for which the state is formed inside the nuclear medium, 
this would indicate that there the interaction depends on its physical size. The 
different fates of the different quarkonium states in cold nuclear matter are due to 
their different sizes, in a hot medium due to the different binding energies. And 
dissociation by color screening should start with a threshold, while the break-up 
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in nuclear matter is presumably continuous. Measurements of charmonia relative 
to open charm in pA up to highest energies (RHIC, LHC) are therefore of great 
importance. 

First applications of the in-medium charmonium study based on the relative survival 
of charmonia vs. open charm were started last year, using LHC data from ALICE and 
CMS [25H28] . In Fig. Ufa), we show mid-rapidity ALICE data for J/ip production at 
intermediate transverse momenta, compared to open charm production in a similar 
kinematic region. In Fig. Hfb), the comparison is extended to larger transverse 
momenta, using CMS data for J/ip production. In both cases, J/ip production 
relative to pp results, scaled by the number of collisions, decreases with increasing 
centrality, as seen by the corresponding R AA values. 
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Figure 4: LHC data from ALICE and CMS [261428] . comparing J/ip production to 
open charm production at intermediate (a) and high (b) transverse momenta 



This decrease has at times been considered as suppressed J / ip production. However, 
that is incorrect: the corresponding R AA for open charm production, as determined 
through D measurements, shows within errors the same behavior. In other words, 
the reduction of the J ftp is in complete agreement with that of open charm; there 
is neither suppression nor enhancement, the fraction of the produced cc pairs going 
into J/ip production has remained in the AA collisions considered here the same as 
in the corresponding pp interactions: 

Raa(J/^) = NAA<yJ l^ = Naa ^ = R AA (cc) (3) 
AA n c N pp (J/ip) n c N pp (cc) AA 

with n c denoting the scaling factor for the number of collisions at the corresponding 
centrality. We therefore have 



N AA {J/j>) = N pp {J/j,) 
N AA (cc) N pp (cc) 



gcc^j/^p- (4) 



ALI-DER-38739 



at all centralities, in the kinematic regime indicated. If we divide the J /if) rates by 
the open charm rates in the same kinematic region (large Pt), this ratio becomes 
centrality independent, and if we normalize it to the corresponding value from pp 
collisions, it becomes unity. The correct and model- independent J/if> survival prob- 
ability for experimental study is thus 



S 



j/il> 



. N AA (cc) . 



/ 



' N pp (j/j>) \ = i ( N AA {j/jjy 

N pp (cc) J g C c^j/^ V n aa(cc) 



(5) 



In the kinematic regime considered so far, it is indeed unity, there is neither sup- 
pression nor enhancement. 

At the LHC, corresponding data for low Pt open charm production is not yet avail- 
able. At RHIC, however, it is provided by the PHENIX and STAR collabora- 
tions [2HHIS], and the relevant comparisons of J/if> vs. open charm production are 
shown in Fig. [5j At high Pt we have a similar behavior as at the LHC, no change of 
J /if) production relative to open charm. At low Pt, however, the R AA (cc) of open 
charm is within errors unity over the whole centrality range; in contrast, R AA (J/if>) 
decreases strongly and thus here gives the correct J /if) survival probability. We see 
that now with increasing centrality, a smaller and smaller fraction of cc pairs go into 
J /if) production, with a suppression of up to 75% for the most central collisions. 
The J /if) production finally surviving could conceivably be largely due to corona 
interactions 
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Figure 5: RHIC data from PHENIX and STAR [30H33] , comparing J /if) production 
to open charm production at high (a) and low (b) transverse momenta 



In this context it seems of interest to note that a decrease of the open charm pro- 
duction rates between central and forward rapidity is observed in both pA and AA 
collisions relative to pp collisions, at Fermilab (pA) and at RHIC (Cu-Cu). We note 
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in particular in the Fermilab data based on 800 GeV pA collisions (see Fig. E]) that 
at forward rapidity there is a suppression of open charm production, relative to 
scaled pp rates [35 J, while at midrapidity there is not. A similar effect is observed 
in Cu — Cu collisions with y/s = 200 GeV at RHIC [36]. If this is passed on to 
charmonium production, one expects for RHIC and LHC data a smaller Raa{J 7V0 
at forward than at mid-rapidity, simply because there are fewer cc pairs there to 
form charmonia, no matter how. Such an effect was indeed observed at RHIC and 
has often been considered a puzzle [37] . 
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Figure 6: The rapidity dependence of open charm production in pA collisions at 800 
GeV, parametrized in the form o v a = A a a pp [35J. 

The aim of this work was to show how quarkonium production has to be calibrated 
in order to obtain experiment-based, model-independent information on suppression 
or enhancement due to the hot medium formed in nuclear collisions. The data 
cited are given for illustration, to make the point. In a more complete study, the 
transverse momentum range of J /if) and open charm should be correlated more 
precisely. The crucial question now remains, of course, how the open charm rates 
for low Pt production at the LHC will behave. That kinematic region is responsible 
for the bulk of cc production, and the issue to be decided is if nuclear collisions do 
or do not lead to an overall reduction. At RHIC, that is not the Fig. Hb) 

shows. 
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